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i i We report on a chiral effective field theory calculation of Compton scattering from the proton. 

Our calculation includes pions, nucleons, and the Delta(1232) as explicit degrees of freedom. It 
uses the "delta expansion", and so implements the hierarchy of scales M K < A = Ma — M N < A x , 
■^j- through a counting M K /A ~ A/A z ~ 8. In this expansion the power counting in the vicinity of 

the Delta peak changes, and resummation of the loop graphs associated with the Delta width is 
indicated. This is designed to extend the region of ^PT applicability to momenta p ~ A. 
We have computed the nucleon Compton amplitude in the delta expansion up to N 3 LO — 
&(e 2 8 A ) — for photon energies (0 ~ M K . This is the first order at which the proton Compton 
scattering amplitudes receive contributions from contact operators which encode contributions to 
the spin-independent polarisabilities from states with energies ~ A x . We fit the coefficients of 
these two operators to the experimental proton Compton data that has been taken in the relevant 
photon-energy domain, and are in a position to extract new results for a p and j3 p . 
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1. Introduction to Compton scattering and Chiral Dynamics 

Chiral dynamics in the baryonic sector is typically thought of as the study of the interactions 
of pions and nucleons. However the dictates of electromagnetic gauge invariance mean that chiral 
symmetry also strongly constrains the interactions of both with photons, and so Compton scattering 
from the nucleon is as fundamental a probe of chiral dynamics as pion-nucleon or nucleon-nucleon 
scattering. The lowest-order term in the Compton scattering amplitude (the long-wavelength limit) 
is the Thomson term which is reproduced by #PT but which, depending as it does only on the 
nucleon charge and mass, is independent of chiral dynamics. However, at shorter wavelengths 
the probing proton starts to be sensitive to the structure of the target. At NLO in heavy baryon 
chiral perturbation theory (HB^PT) the dominant new contribution comes from a single pion loop 
with photons coupling to the pion or to the 7rN vertex (see Fig. 1), and hence a prediction can be 
made for these structure effects. This includes, but is not limited to, the numbers known as the 
polarisabilities of the nucleon; the latter are the first terms in an expansion of the six scattering 
amplitudes in powers of the photon energy. 

The application of chiral dynamics to Compton scattering dates right back to the dawn of 
baryon #PT, and most famously the lowest-order predictions for the electric and magnetic polaris- 
ability of the proton in HB^PT, a = 12.5 x 10~ 4 fm 3 and jS = 1.2 x 10~ 4 fm 3 [1] are in extremely 
good agreement with experimental determinations. Higher-order calculations [2, 3], and calcula- 
tions of spin polarisabilities [4, 5, 6, 7, 8] followed. The first studies to systematically compare the 
full fourth-order predicted cross-section to a compendium of experimental data were published by 
some of the current authors [9, 10], first with values of a and /3 taken from the Particle Data Group, 
then using them as fit parameters (at fourth order, they are no longer predicted but have contribu- 
tions from LECs). These studies obtained an excellent fit at low energies, and clearly demonstrated 
the pion-production cusp. Furthermore the fourth-order results are in markedly better agreement 
with the data than the third-order ones (even when the freedom to fit the polarisabilities is not ex- 
ploited). However they also showed that the cross-section completely failed to keep pace with the 
sharp rise in the experimental data above the cusp, particularly at backward angles. This is not 
surprising: from around 200 MeV the experimental cross section is dominated by the Delta reso- 
nance. While contributions from the Delta enter the Delta-less theory via the LECs, these clearly 
cannot describe the resonance. Nonetheless, by restricting the fit to appropriate regions of energy 
and angle, the following results were obtained [10]: 

a p = (12.1 ± 1.1 (stat.))+^ (theory) x 10~ 4 fm 3 
f3 p = (3.4 ± 1.1 (stat.))±8:i (theory) x 10~ 4 fm 3 . 

Work has been done on including the Delta in chiral effective theories from the earliest days 
of baryon #PT, motivated partly by its role as the large-Afc degenerate partner of the nucleon, and 
partly by its experimental relevance. However there is no universal agreement as to the correct 
power-counting; it is not obvious how to treat the scale A = Ma — M N ps 2M n . The most commonly 
adopted solution is to treat them as proportional to the same scale (the "small scale expansion"), 
a pragmatic treatment which however is not strictly in the spirit of #PT as A does not vanish 
in the m q — > limit. With this treatment 7iA loops come in at the same order as nN loops. In the 
particular case of Compton scattering the lowest-order prediction for a and j8 is destroyed with this 
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counting. More recently Pascalutsa and Phillips proposed a more sophisticated treatment in which 
the power counting differs in different energy regions, depending on which scales are enhanced. 

At low energies they introduce the ratio 8, with 5 = ~ x~ so that <5 2 = ( ^- J ~ j^. With this 
"delta-expansion", the nA loops enter at one order higher than nN loops. 

Compton scattering cross sections calculated at lowest non-trivial order in the small scale 
expansion were compared to data by Hildebrandt etal. [11]. Because of the very large contribution 
of the Delta to a and j3 , they had to introduce counterterms for these, which they then fit to data 
in a similar fashion to Ref. [10]. Interestingly, in spite of the rather different input and data set, the 
results were very similar: 

a p = 11.5±1.4 stat x l(T 4 fm 3 , p p = 3.4 ± 1.6 stat x l(T 4 fm 3 . 

However the Delta-full results were able to describe data up to around 200 MeV lab energy, unlike 
the Delta-less results. 

Clearly the next step is to include both the Delta and higher-order effects, and that is the aim 
of the current program. Very recently Lensky and Pascalutsa have also taken steps in this direction, 
using a relativistic formulation. However they have not reached the order at which they must fit a 
and j8, and they have so far only compared to a subset of world data [12]. 

2. Comments on the proton Compton database 

Modern experiments measure the differential cross section for tagged photons over a range 
of energies and angles. Between 1956 and 1995, twelve Compton scattering experiments took 
data at energies below 200 MeV. In a useful paper Baranov et al [13] examined the world data on 
proton Compton scattering, and demonstrated that the data from the 50s and 60s was compatible 
with the more modern data from 1974 onwards, and was useful in reducing errors. In addition, the 
most recent and extensive experiment used the TAPS detector and MAMI tagged photon facility at 
Mainz [14]. 

As part of the current project we have evaluated the consistency of the various data sets in as 
model-independent a fashion as we can. We find that it is crucial to allow the normalizations to 
float (as done in Ref. [13, 10] but not in [11]). For the Mainz data we increase the statistical error 
by including a "random systematic" error of 5% in quadrature, as described in Ref. [14]. Some of 
the very early data cannot be accommodated, nor — as is well known — can the Bonn data of Genzel 
et al. [15]. With these caveats, we find that the world data is reasonably consistent in the range 
0-240 MeV. There is enough noise in the data sets that a % 2 of one per degree of freedom is not 
obtainable though. 

Such data has traditionally been analyzed in one of two ways. Early experiments typically 
took very low energy data (up to around 100 MeV) and fitted to a low-energy expansion of the 
cross section which included only Born contributions and the electric and magnetic polarisabilities, 
with these as fit parameters. However at these energies the Born cross section dominates, and 
the sensitivity to a and j8 is limited. A desire to use the full range of data up to the region of 
the pion-production threshold meant that later authors turned to dispersion relations to obtain a 
prediction for the variation of the cross section with energy, leaving only certain combinations of 
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Figure 1: Sample diagrams that contribute at &(8 A ). 

the polarisabilities as free parameters to be fit. Olmos de Leon et al. [14], for instance, fit a, /3 and 
Yjc. Unfortunately the dispersion relation approach is not entirely free of assumptions; for instance 
the asymptotic form of the backward spin-independent amplitude is assumed to be dominated by 
a exchange. There is only one attempt, by Baranov et al. [13] to fit the whole data set (up to 
150 MeV) which was unfortunately done before the Mainz data was released. 

The use of HB^PT with or without the Delta provides an alternative, model-independent ap- 
proach to the problem, and the history of such attempts is given above. 



3. An analysis of proton Compton data using the delta-expansion 

In figure 1 we show the type of diagrams which contribute at each order. Most of the param- 
eters which enter are well known. The values of the second order 7iN scattering LECs c, which 
we use must of course be ones appropriate for a Delta-full theory. The yNA coupling constant 
(variously b\ or Gm in the literature) is not well known; though it has been fit to resonant elec- 
troproduction, the power counting in that regime differs from that appropriate to the low-energy 
regime, and so it is not clear that the same value should be used. We take the same approach as 
Hildebrandt et al. [11] and include it as a fit parameter. 

There is a problem with the fit when including both 7iA loops and higher-order TiN terms. Both 
individually correct the tendency of the third-order kN cross section to be too low in the vicinity 
of the cusp at backward angles. Together they over-correct. It is easy to deduce that the reason is 
related to the spin polarisabilities: there is exactly the same over-correction in the value of Jmi- We 
therefore promote a counterterm from higher order to allow us to set this to the best estimates of its 
experimental value [16, 11] or to fit it to data. 

There is an obvious tension in deciding how high to take our energy cut-off in fitting data. 
On the one hand the higher we take it, the more data we can include and the lower the statistical 
error. On the other hand we know that we are not using the correct power counting for the Delta 
resonance region and we would not expect a good description of the data there. There are also 
issues associated with the choice of frame for the calculation which become more acute at higher 
energies. We have looked at cut-offs up to 240 MeV lab energy. 

At this point our fits are preliminary, and we will not quote results. Figure 2 however shows 
a fit demonstrating that we can obtain a good description of data up to at least 200 MeV. The two 
curves illustrate the effect of differences in the treatment of the Delta: in one the treatment strictly 
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Figure 2: Fits to world data: with (green, lower) and without (red, higher) the inclusion of the Delta width 
and higher order terms in the yNA vertex. The data are from a wide variety of sources which we will not list 
here. For display purposes they have been binned in 40° bins. 

follows the low-energy power counting, in the other we have included the Delta width and A/M N 
corrections to the yNA vertex. The parameters used are consistent with those obtained by Beane et 
al. [10] and Hildebrandt et al. [1 1]. In particular a largish value of /3 is preferred, a result consistent 
with the recent analysis of Lensky and Pascalutsa [12] 

A main goal of this study is to produce reliable single-nucleon amplitudes for incorporation in 
a code for scattering from the deuteron [17], which will allow the analysis the data that is currently 
being taken at MAX-lab [18]. We believe we have essentially reached that goal. However the 
prospects for reducing the error bars on a p and fi p are not so good. More high-quality lower-energy 
data (below pion production) is urgently needed. 
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